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ATTENUATION OF SINUSOIDAL PERTURBATIONS SUPERIMPOSED 


ON LAMINAR FLOW OF A LIQUID IN A LONG LINE 
by Carl M. Holland, Robert J. Blade, and Robert G. Dorsch 
Lewis Research Center 


SUMMARY 

The attenuation constant for sinusoidal pressure and flow perturbations 
superimposed on the laminar flow of a viscous liquid was measured in a system 
consisting of a long, straight, cylindrical hydraulic line. The upstream and 
downstream ends of the line were securely fastened to the ground. A sinusoidal 
perturbation was imposed on the mean flow at the upstream end by means of a 
small oscillation of a throttle valve about a partly open mean position. The 
downstream end was terminated in a restricting orifice. Pressure perturbations 
were measured at three locations along the line for frequencies from 15 to 
100 cps. These pressure measurements were reduced by use of a pair of complex 
damped acoustic one -dimensional wave equations to obtain the attenuation con- 
stant along with the phase constant and the dimensionless downstream admittance. 
For the range of frequencies investigated, the experimental values of the atten- 
uation constant are in good agreement with classical theory. 


INTRODUCTION 

Acoustic pressure and flow perturbations in hydraulic lines are of impor- 
tance in the operation of various fluid systems. The dynamics of hydraulic con- 
trol and rocket -propellant feed systems are of particular interest to the aero- 
space industry. For example, longitudinal instability of liquid rocket boosters 
results from dynamic coupling between flow perturbations in the propellant feed 
system and structural oscillations. The Lewis Research Center has therefore 
conducted experimental and analytical studies of acoustic disturbances in 
flowing- liquid systems. 

Transfer functions for acoustic disturbances superimposed on the turbulent 
flow of a liquid propellant flowing in long cylindrical lines of several geomet- 
rical networks are reported in references 1 to 3. These measurements were made 
with a nonviscous hydrocarbon liquid with small attenuation (about 3 percent 
per 100 ft). Other investigators have reported transfer functions and distrib- 
uted line impedance for damped acoustic disturbances in lines containing a vis- 
cous liquid with appreciable attenuation (refs. 4 to 7). 

The analytical expression for the pressure transfer function depends upon 



three parameters: the attenuation constant, the phase constant, and the ratio 

of the downstream admittance to the characteristic admittance of the line 
(referred to herein as the "dimensionless downstream admittance"). Experimental 
values of the attenuation constant are reported in references 8 and 9 for sound 
waves in stationary air within tubes. 

The results of an experimental study of the attenuation constant for acous- 
tic disturbances superimposed on the mean laminar flow of a viscous liquid flow- 
ing in a long line are reported herein. Experimental data were obtained by 
measuring perturbation pressures at three stations along the line. Experimental 
values were computed from the data and compared with analytical values obtained 
from classical theory. 


APPARATUS 

Flow system. - The essential parts of the open- loop pumped-return flow 
system used in the experiment are shown in figure 1. The test fluid was an 
additive-free SAE 20 motor oil with a viscosity which resulted in laminar mean 
flow in the line at the operating conditions. The oil was supplied by a gear 
pump, and the mean flow rate was measured by means of a rotameter. A heat ex- 
changer in the fluid-supply tank capable of either heating or cooling maintained 
the fluid at a constant temperature during each run. Hydraulic accumulators 
were placed between the pump and the test line to provide steady supply pressure. 
The discharge from the test line was submerged in a constant -height, vented tank. 
Fluid was returned to the supply tank by intermittent operation of the return 
pump. 


Flow disturbance generator. - Sinusoidal perturbations of flow and pressure 
were induced in the system by means of an electrohydraulic servoactuated throt- 
tle located just upstream of the test line. The throttle was oscillated sinus- 
oidally about a partly open mean position in response to an alternating voltage. 

Test line. - The test line, made of stainless-steel tubing, 1.00 inch in 
outside diameter with a 0.065-inch wall thickness, was 68 feet long. A sche- 
matic diagram of the test line is shown in figure 2. The upstream end was 
attached to the throttle valve, which, in turn, was securely fastened to the 
ground. The downstream end was bolted to a 3- inch-diameter pipe fastened to a 
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5 -ton block of concrete set in the ground. Between the ends, the test line 
rested on horizontal transverse wires . spaced at 2-foot intervals. The line was 
terminated in an orifice plate containing 21 holes , 0.040 inch in diameter. 

The orifice plate was rigidly attached to the line. 

Instrumentation . - Pressure sensors, commercial flush diaphragm units , 
were rigidly attached to the line at stations A, B, and C (fig. 2). The pres- 
sure sensors with their associated amplifiers were calibrated by static methods. 
Relative gain factors among the sensors were then determined from the calibra- 
tions. 

The output of the pressure sensors was in the form of alternating-current 
electrical signals with phase and amplitude determined by the local perturbation 
pressure. These alternating-current signals were amplified and fed into a 
transfer -function analyzer. The transfer -function analyzer consisted of two 
units: (l) an oscillator, which supplied (a) the signal to activate the electro- 

hydraulic flow-disturbance generator and (b) a four -phase reference voltage to 
the resolved component indicator, and (2) a resolved component indicator, which 
indicated the in-phase and quadrature (90° out of phase) components of the pres- 
sure signals with respect to the reference voltage. The transfer-function 
analyzer effectively rejects all frequencies (associated with noise and har- 
monics) except the reference frequency (generated by the oscillator) and indi- 
cates the output signal of the pressure sensors in resolved component form. 


PROCEDURE 

Test. - The flow system was operated at constant mean flow rate, pressure, 
and temperature until conditions had stabilized. Mean gage pressure readings 
were taken at stations A and C, the mean flow rate was read on the rotameter, 
and the temperature adjacent to the downstream orifice was recorded. The servo- 
throttle was operated at a series of frequencies. The amplitude of the throttle 
area variation was maintained constant over the entire range of frequencies for 
each run. It was kept small relative to the mean open area to avoid nonlinear 
effects but large enough to give an adequate ratio of signal to noise. The 
resulting average amplitude of sinusoidal pressure perturbations was about 
15 pounds per square inch. At each frequency, the values of the in-phase and 
quadrature components of the three dynamic pressures (upstream, quarter length, 
and downstream) were read in succession on the resolved component indicator. 
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The mean pressure drop across the downstream orifice and the fluid temper- 
ature were constant during each run hut were independently changed from run to 
run to vary the average orifice impedance (controlled by both of these factors) 
and the viscosity of the fluid (controlled mainly by fluid temperature). 

The data were taken over a range of disturbance frequencies from 15 to 
100 cps, mean line pressures from 165 to 315 pounds per square inch gage, and 
temperatures from 80° to 125° F. The range of variation of these parameters de- 
pended upon facility limitations. The range of frequencies, mean line pres- 
sures, and temperatures resulted in a range of absolute viscosities from 
0. 748x10" ^ to 2. 5 4x10" ^ pound-second per square foot, Reynolds numbers of the 
mean flow (with inner tube diameter as the characteristic length) from 263 to 
1200, mean flow rates from 5.5 to 8.0 feet per second, and magnitudes of dimen- 
sionless downstream admittance from 0.63 to 1.2. 

The oil was deaerated between runs by flowing it at low speed and high 
temperature for about 1 hour through the receiving tank (fig. 1, p. 2) in which 
a vacuum of 10 inches of mercury existed. During each data run, the receiving 
tank was vented to the atmosphere. 

Data analysis . - The following quantities were computed from the experi- 
mental data for each frequency: 

(1) The attenuation constant a in nepers per foot 

(2) The phase constant |3 in radians per foot 

(3) The complex dimensionless admittance at station C, Z Q / Zq 
(A ll symbols are defined in appendix A. ) 

These quantities were computed from the data for each frequency by use of 
a pair of complex pressure transfer functions for (l) stations A and C (fig. 2) 
and (2) stations B and C. This pair of transfer functions is given by equa- 
tions (l) and (2), respectively: 

Pa , N Z Q / 

— = cosh(a + jp)Z A + — sinh(a + j(3)Z A (l) 


- cosh(a + jp)Z-D + sinh(a + j p ) Z 


B 


B 


( 2 ) 


These equations were obtained from the well-known complex pressure transfer 
function for damped sinusoidal one -dimensional acoustic waves in uniform pipes. 
They were solved (see appendix B) on a digital computer by use of a numerical 
method. 


The input data necessary were (l) the measured values of the in-phase and 
quadrature components of P A , Pg, and Pq and (2) the distances in feet Za 
and Zg between stations A and C and between stations B and C, respectively. 
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square foot (0.0358 N-sec/sq m); Reynolds number of mean flow, 1200. 



square foot (0.070 N-sec/sq m); Reynolds number of mean flow, 431. 



(c) Run c. Absolute viscosity of fluid, 2. 54x10"^ pound-second per 
square foot (0.122 N-sec/sq m); Reynolds number of mean flow, 263. 


Figure 3. - Variation of attenuation constant with frequency. 
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For comparison, theoretical values of the attenuation constant a, the 
phase constant [3, and the specific characteristic impedance of the line Z Q 
were computed from classical equations derived in appendix C (primarily equa- 
tions (C16), (C17), and (C18), respectively). The values of the density and 
the effective adiabatic hulk modulus of the liquid in the line required by these 
equations were obtained from the manufacturer and from references. (See 
appendix D for these values and references). The absolute viscosity in pound- 
seconds per square foot was obtained for each run from the mean-state data 
(mean pressures at stations A and C and mean flow rate) by use of Fbiseuille f s 
coefficient (ref. 10, p. 249). 


RESULTS AND DISCUSSION 


Experimental and analytical results for three runs (denoted herein as a, 
b, and c) that cover the range of mean line conditions are given in figures 3 
to 5. The data points were calculated from experimental pressure measurements 
by the methods described in the PROCEDURE section and the analytical curves were 
obtained from the equations of appendix C. 



(a) Run a. Absolute viscosity of fluid, 0. 748xl0 -3 pound-second per 
square foot (0.0358 N-sec/sq m); Reynolds number of mean flow, 1200. 



(b) Run b. Absolute viscosity of fluid, L 46x10“ 3 pound-second per 
square foot (0. 070 N-sec/sq m); Reynolds number of mean flow, 431. 



(c) Run c. Absolute viscosity of fluid, Z 54x10 3 pound-second per 
square foot (0. 122 N-sec/sq m); Reynolds number of mean flow, 263. 

Figure 4. - Effect of frequency on attenuation, expressed as pressure 
ratio e _100a , where a is attenuation constant 


Figure 3 gives the atten- 
uation constant a in nepers 
per foot as a function of fre- 
quency. Dashed data-fit curves 
are shown along with the analyt- 
ical curves and the data points. 
The data-fit curves were ob- 
tained by a least-squares 
method (appendix E) in which the 
attenuation constant was taken 
to be proportional to the square 
root of frequency. This func- 
tional relation was chosen for 
the data-fit curves because it 
is commonly used as an approxi- 
mation (ref. 10, p. 242) to the 
exact analytical curve. Al- 
though there is considerable 
scatter in the individual data 
points, the data-fit and analyt- 
ical curves are considered to 
agree within the limits of ac- 
curacy of the experiment. 

As an illustration of the 
effect of the attenuation con- 
stant on acoustic disturbances 
traveling down the line, fig- 
ure 4 presents the quantity 
e -100a as a function of fre- 
quency. This quantity is the 
ratio of the amplitude of the 
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(a) Run a. Absolute viscosity of fluid, 0. 748x10’^ pound-second per 
square foot (0.0358 N-sec/sq m); Reynolds number of mean flow, 1200. 
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(b) Run b. Absolute viscosity of fluid, L 46x10 pound-second per 
square foot (0.070 N-sec/sq m); Reynolds number of mean flow, 431. 
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(c) Run c. Absolute viscosity of fluid, 2. 54x10“^ pound-second per 
square foot (0. 122 N-sec/sq m); Reynolds number of mean flow, 263. 

Figure 5. - Effect of frequency on ratio of phase constant to frequency. 


downstream to the upstream pressure perturbation for a wave traveling in the 
downstream direction for two points along a line which are separated by 100 feet 
(30.5 m). Figure 4 shows that for the conditions of this experiment the atten- 
uation of pressure amplitude was of the order of 25 percent in a 100-foot line. 

In addition to the attenuation constant , experimental and analytical values 
of the phase constant (3 (see eqs. (l) and (2)) were determined. These values 
are presented in figure 5 in the form p/f as a function of frequency. The 
experimental and analytical results are in good agreement. The quantity p/f 
is inversely proportional to the velocity of propagation of sinusoidal disturb- 
ance waves in the line c. More precisely, 

f-T < 3) 

Figure 5 indicates that the value of the velocity of propagation was nearly 
constant for most frequencies of measurement. 

As a check against obtaining extraneous roots in computing a, the experi- 
mental magnitude and angle of the complex dimensionless admittance at station C 
were compared with semianalytical values (see appendix 33). These experimental 
and calculated values corresponding to the runs shown in figures 3 to 5 are 
given in table I. The mean experimental values of the magnitude of this ratio 
were in very good agreement with the calculated values; the maximum deviation of 
experiment from the calculations was less than ±11 percent. The maximum differ- 
ence between the experimental and the calculated angles of this ratio was less 
than ±10°. This agreement confirms the implied assumption that the method of 
computing the data gave the roots of interest. 
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TABLE I. - CALCULATED AND EXPERIMENTAL VALUES OF COMPLEX DOWNSTREAM DIMENSIONLESS ADMITTANCE Z o /Z 0 



Run 

a 



Run b 

— 

Run c 

Frequency, 

cps 

Magnitude of 
admittance 
2 0 / z C 

Angle of 
admittance 
z o/ z C-> 
deg 

Frequency, 

cps 

Magnitude of 
admittance 

Zo/Zc 

Angle of 
admittance 
z o/ z C-> 
deg 

Frequency, 

cps 

Magnitude of 
admittance 

z 0 /z c 

Angle of 
admittance 

z o/ z C' 

deg 

Experi- 

mental 

Calcu- 

lated 

Experi- 

mental 

Calcu- 

lated 

Experi- 

mental 

Calcu- 

lated 

Experi- 

mental 

Calcu- 

lated 

Experi- 

mental 

Calcu- 

lated 

Experi- 

mental 

Calcu- 

lated 

17 

0.65 

0.67 

0 

-2 

12 

1.11 

1.08 

-7 

-4 

22 

1.10 

1.07 

-2 

-4 

22 

.65 

.67 

0 

-2 

14 

1.08 

1.08 

-5 

-4 

24 

1.08 

1.07 

-3 

-4 

24 

.65 

.67 

0 

-2 

17 

1.12 

1.07 

-6 

-3 

27 

1.06 

1.07 

-1 

-3 

27 

.64 

.66 

0 

-2 

22 

1.12 

1.07 

-4 

-3 

29 

1.04 

1.07, 

-1 

-3 

29 

.64 

.66 

0 

-2 

24 

1.12 

1.06 

-2 

-3 

32 

1.01 

1.06 

-1 

-3 

32 

.63 

.66 

0 

-2 

27 

1.07 

1.06 

-1 

-3 

34 

1.00 

1.06 

-2 

-3 

34 

.63 

.66 

1 

-2 

29 

1.05 

1.06 

0 

-3 

37 

.95 

1.06 

-4 

-3 

39 

.64 

.66 

0 

-2 

32 

1.00 

1.06 

2 

-2 

39 

.96 

1.06 

-6 

-3 

42 

.64 

.66 

-4 

-1 

34 

1.01 

1.06 

0 

-2 

42 

1.02 

1.06 

-10 

-3 

44 

.66 

.66 

-5 

-1 

37 

1.02 

1.05 

-1 

-2 

44 

1.03 

1.06 

-8 

-3 

47 

.67 

.66 

-3 

-1 

39 

.95 

1.05 

-4 

-2 

47 

1.08 

1.05 

-8 

-3 

52 

.68 

.66 

1 

-1 

42 

1.01 

1.05 

-10 

-2 

49 

1.09 

1.05 

-6 

-3 

54 

.69 

.66 

1 

-1 

44 

.99 

1.05 

-9 

-2 

52 

1.10 

1.05 

-5 

-3 

57 

.67 

.66 

3 

-1 

47 

1.05 

1.05 

-8 

-2 

54 

1.10 

1.05 

-5 

-3 

62 

.66 

.66 

3 

-1 

49 

1.13 

1.05 

-3 

-2 

57 

1.07 

1.05 

-4 

-2 

67 

.66 

.66 

4 

-1 

52 

1.10 

1.05 

-4 

-2 

59 

1.07 

1.05 

-4 

-2 

69 

.63 

.66 

4 

-1 

57 

1.10 

1.05 

-3 

-2 

62 

1.02 

1.05 

-3 

-2 

72 

.63 

.66 

0 

-1 

62 

.99 

1.04 

1 

-2 

64 

1.05 

1.05 

-4 

-2 

74 

.64 

.66 

1 

-1 

64 

1.03 

1.04 

-3 

-2 

67 

.99 

1.05 

-6 

-2 

77 

.64 

.66 

-2 

-1 

67 

1.02 

1.04 

-4 

-2 

69 

1.00 

1.05 

-5 

-2 

79 

.68 

.66 

-4 

-1 

69 

1.00 

1.04 

-6 

-2 

72 

1.00 

1.05 

-6 

-2 

82 

.72 

.66 

-2 

-1 

72 

1.03 

1.04 

-11 

-2 

74 

1.00 

1.05 

-8 

-2 

84 

.70 

.66 

1 

-1 

74 

1.06 

1.04 

-9 

-2 

77 

1.04 

1.04 

-7 

-2 

87 

.71 

.65 

3 

-1 

77 

1.05 

1.04 

-7 

-2 

79 

1.11 

1.04 

-8 

-2 

89 

.70 

.65 

7 

-1 

79 

1.06 

1.04 

-7 

-2 

82 

1.11 

1.04 

-5 

-2 

92 

.69 

.65 

8 

-1 

82 

1.08 

1.04 

-5 

-2 

84 

1.04 

1.04 

-5 

-2 

94 

.67 

.65 

7 

-1 

84 

1.09 

1.04 

-5 

-2 

87 

1.06 

1.04 

-4 

-2 

97 

.68 

.65 

-4 

-1 

87 

1.08 

1.04 

-4 

-1 

89 

1.06 

1.04 

-5 

-2 






89 

1.06 

1.04 

-4 

-1 

92 

1.06 

1.04 

-6 

-2 






92 

1.05 

1.04 

-6 

-1 

97 

1.09 

1.04 

-7 

-2 






94 

1.03 

1.04 

-4 

-1 

99 

1.09 

1.04 

-7 

-2 


97 

99 

102 


.99 

.97 

.95 


1.04 

1.04 

1.04 


-3 

-6 

-7 



CONCLUSIONS 


In a study of the attenuation of sinusoidal perturbations superimposed on 
laminar flow of a liquid in a long line, the experimental and analytical results 
were in good overall agreement. These results indicate that one -dimensional 
damped acoustic -wave theory, as presented herein, satisfactorily describes the 
attenuation constant and, more generally, sinusoidal disturbance propagation. 
This statement is restricted to laminar mean flow of viscous Newtonian liquids 
within long, firm- walled, nonvibrating, straight, cylindrical pipelines within 
the frequency range studied. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, September 9, 1965. 
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APPENDIX A 


SYMBOLS 


A constant defined ty eq. (E4) 
a tube inner radius , ft (m) 

c velocity of sinusoidal wave propagation in line, ft/sec (m/sec) 
f frequency, cps (hertz) 

I defined by eq. (CIO) 

J q Bessel function of first kind of order zero 

J^_ Bessel function of first kind of order one 

j -vA 

K effective adiabatic bulk modulus of liquid in line, lb/sq ft (w/m 2 ) 
l distance from station C (fig. 2), ft (in) 

W defined "by eq. (C9) 

n number of data points in run 

P sinusoidal-disturbance complex pressure amplitude, P(x), lb/sq ft (h/m 2 ) 

P q sinusoidal-disturbance complex pressure amplitude at x = 0 (station 0) 

for component of wave that is traveling in negative x-direction, 
lb/sq ft (h/m 2 ) 

p disturbance pressure, p(x,t), lb/sq ft (h/m 2 ) 

R specific series resistance of line to sinusoidal perturbation, 
lb-sec/cu ft (N-sec/m 3 ) 

Re Reynolds number of mean flow with use of inner tube diameter as 
characteristic length 

r radial coordinate in any cross-sectional plane of line with r = 0 
along axis, ft (m) 

t time, sec 

U complex amplitude of sinusoidal-perturbation velocity in x-direction 
averaged over cross-sectional area of line, U(x), ft/sec (m/sec) 
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F 


U’(x,r) complex amplitude of sinusoidal-perturbation velocity in x-direction, 
ft/sec (m/sec) 

u perturbation velocity in x-direction averaged over cross-sectional 

area of line, u(x,t), ft/sec (m/sec) 

u ! perturbation velocity in x-direction, u ! (x,r,t), ft/sec (m/sec) 

v* perturbation velocity in r-direction, v’(x,r,t), ft/sec (m/sec) 

x coordinate in axial direction of line with x = 0 at station C 

(fig. 2) and positive direction in upstream direction, ft (m) 

Z q specific impedance at station C (fig. 2) looking downstream, a complex 

number (approximately real in this experiment), lb-sec/cu ft 
(N-sec/m^) 

Z G specific characteristic impedance of line, a complex number, 

lb-sec/cu ft (N-sec/m 3 ) 

a attenuation constant, Np/ft (Kp/m) 

(3 phase constant, rad/ft (rad/m) 

{x absolute viscosity of fluid, lb-sec/sq ft (N-sec/m 2 ) 

P density of fluid (denotes mean value except in dp), slugs/cu ft (kg/m 3 ) 

p e effective density of fluid defined by eq. (C12), slugs/cu ft (kg/m^) 

cn frequency, 2 rtf, rad/sec 

Subscripts : 

A station A (fig. 2) 

B station B (fig. 2) 

C station C (fig. 2) 

F data fit 

i data point 

x real part of complex number 

y imaginary part (quadrature component) of complex number 
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APPENDIX B 


DETERMINATION OF PRESSURE TRANSFER FUNCTION PARAMETERS FROM DATA 

Using identities and equating the real and the imaginary components of 
equations (l) and (2) gives four equations in the four unknowns a, (Z q /Zq) x , 

and ( Z g /Zq) as follows: 

y 

= (cosh ctZ A )cos pZ A + (sinh a,Z A )cos pZ A - ' (cosh aZ A )sin pZ A 


(Bl) 


= (sinh aZ A )sin PZ A + (cosh aZ A )sin pZ A + J (sinh aZ A )cos pZ A 

y ° x y 



(B2) 


/P B \ fz 0 \ /z 0 \ 

= (cosh aZ B )cos pZ B + y (sinh aZ B )cos pZ B - (cosh a.Z B )sin pZ B 


(B3) 


= (sinh a,Z B )sin pZ B + (cosh aZ B )sin pZ B + (^~j (sinh aZ B )cos pZ B 


(B4) 


Eliminating (Z d /Zq) and (Z c /Zq) from equations (Bl) to (B4) and using 

x y 

identities gives two equations in the two unknowns a and p : 


(sinh aZ A )cos pZ A 

/ p a\ 

^ 0 / 
X 


(cosh aZ A )sin PZ A j 

i'PaN 

pc/ 

y 

(sinh aZ B )cos pZ B | 

' p b\ 

< p ci 


(cosh aZ B )sin pZg | 

' p b\ 

iP 4 


- [sinh a(Z A - Z B )]cos p( z A - Z B ) = 0 


(B5) 
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r 



(cosh aZ A )sin (3Z A 


© 

(sinh aZ A )cos 



_ 

X C/ y 



(cosh aZ B )sin 0Z B 


(l) 

(sinh aZ^)cos 

P Z B 



c/ y 




-[sin (3 (Z A - Z B )] cosh a(Z A - Z B ) = 0 (B6) 


Equations (B5) and (B6) were solved for a and (3 on a digital computer by 
use of a numerical solution^ the Newton-Raphson method. The criterion of con- 
vergence both a and (3 were forced to satisfy was that the magnitude of (l) 
the difference between two successive values of the unknown (a, (3) computed by 
the numerical method for a given frequency divided by (2) the last computed 
value of the unknown at that frequency must be less than or equal to 5x10“^. 
Initial estimates for a and (3 were made for the calculation at the highest 
frequency of each run. The values to which a and (3 converged at a given 
frequency were used as the initial estimate in the numerical solution at the 
next lower frequency of measurement. 

The unknowns (Z 0 /Z n) and (Z n /Zn) can be calculated as follows: Equa- 

tions (Bl) and (B2) give X y 


= (pq) ( sinh aZ A^ cos “ (p^) ( GOSh a ^A^ sln 


■][• 


+ (cos |3Z A )sin pZ A | |( cosh' 2 aZ A ) sin 2 0Z A + ( sinh 2 aZ A ) cos'^pz^ 


a 


-l 


and 


(B7) 


( 1 ) -[( 

Tr ) ( sin P^a) cos1:i aZ A + ( t 

^ ) (cos pZ A ) sinh a i A 

\ / x L\ 

c/ y V 

c 4 
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- (sinh aZ A )cosh a Z A I I ( cosh 2 aZ A ) sin 2 (3Z A + ( sinh 2 aZ A ) cos 2 (3Z A 


' -1 

A 


(B8) 


The values of (Z 0 /Zq)_^ and (Z q /Zq)_^ were computed for each frequency by sub- 
stituting the values of a and (3, obtained from the numerical solution of 
equations (B5) and (B6)^ into equations (B7) and (B8). Note that the experi- 
mental values of (P A /Pc) x and - ( p a/ p c)^ are explicitly required in equa- 
tions (B7) and (B8). 
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As a check against obtaining extraneous roots in computing a, the experi- 
mental magnitude and angle of the complex dimensionless admittance at station C 
were compared with values obtained from (l) an equation for the specific char- 
acteristic impedance of the line (eq. (C18)) and (2) an approximation of the 
specific impedance at station C looking downstream (fig. 2) obtained by taking 
it as the resistance of the downstream orifice found from the slope of a curve 
of mean pressure against fluid velocity. These experimental and calculated 
values corresponding to the runs shown in figures 3 to 5 are given in table I 

(p. 8). 



APPENDIX c 


DERIVATION OF ANALYTICAL PRESSURE TRANSFER FUNCTION AND RELATED PARAMETERS FOR 
SMALL SINUSOIDAL PERTURBATIONS SUPERIMPOSED ON MEAN LAMINAR FLOW 

Assumptions 

The following assumptions lead to a system of partial differential equa- 
tions : 

(1) The fluid was considered Newtonian. 

(2) Stokes hypothesis (ref. 11, p. 383) was assumed to apply. 

(3) The medium was continuous and uniform. 

(4) The absolute viscosity was considered constant. 

(5) Axial symmetry was assumed to exist, because the line was cylindrical 
and its axis was straight. 

(6) The perturbation flow was regarded as laminar. 

(7) The viscous force that is proportional to S 2 u'/Sx^ (ref. 5) and end 
effects were neglected, because the ratio of tube inner radius to line length 
was small, 6.1xl0“ 4 . 

(8) The pressure across the cross section of the fluid can be considered 
constant, because of (a) the existing condition j_i/a-s/pK « 1 (ref. 7), 

(b) the rigidity of the pipe walls (ref. 6), and (c) the straightness of the 
pipe. 


(9) The nonlinear convective acceleration terms in the momentum equation 
were neglected, because the ratio of mean flow speed to the velocity of 
sinusoidal-disturbance wave propagation was small, less than 5xl0“ 3 (ref. 5). 

(10) The change of density with position at any instant can be neglected, 
because the adiabatic bulk modulus of the liquid in the line was at least 
3.0x10^ pounds per square foot (ref. 5). 

(11) The energy equation had little effect on the solution, because (a) 
heat conduction can be ignored because of the small difference between the 
ratio of specific heats and 1 (0.159) and the large Rrandtl number (about 28) 
(refs. 10, p. 243, and 7) and (b) heat radiation can be ignored. 

(12) The classical thermodynamic equation of state was used, because atten- 
uation due to molecular exchanges of energy (discussed in ref. ll) was assumed 
negligible. 
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Basic Equations 


It follows from the assumptions that the system is described by three 
partial differential equations (ref. 5). The momentum equation is 


dp du * ll d 

53E = p 5T" ' r dF 



(ci) 


The continuity equation is 


dp 

dt 


+ 


dv' 
dr" '' 




du' 

dx 


= 0 


(C2) 


The thermodynamic equation of state for a liquid is 


dp dp 
P ” K 


(C3) 


Combining equations (C2) and (C3) and averaging the result over the cross- 
sectional area of the line gives (ref. 5) 




(C4) 


The steady- state solution only is sought. For a sinusoidal disturbance, the 
perturbation velocity in the axial direction u* and the perturbation pressure 
p can be written as 


u* = U ! (x,r )e (C5) 

p = P eJart (C6) 

where P is a function of x only. For a tube wall that is stationary in the 
axial direction, the component of fluid velocity in the axial direction at the 
wall of the tube is zero. Therefore, 


U f (x,a) = 0 


(C7) 


Substituting equations (C5) and (C6) into equation (Cl), solving for U ! (x,r) 
with the use of equation (C7), and averaging the resulting equation over a cross- 
sectional area gives the complex amplitude of the velocity perturbation averaged 
over a cross-sectional area as (ref. 10, p. 240) 


where 


U = 


dP 

3x_ 

l_iN 2 


£ 


2 Ji(aN)\ 

aN J 0 (aW ) ) 


H = (1 



(C8) 


(C9) 


16 



f 


Define 


Rearranging equation (C8) gives 


_ 2 JlM 

a J 0 (aW) 


“ = ( R + J“P e ) U 


where the effective density of the fluid (also referred to as the 
specific series inertance" ) is 



and the distributed specific series resistance is 

- PoVfp Ux + x y) 


R = 


(vf ■ ■■) *(«•'.) 

Generally following the method of reference 10 (pp. 233-236), that 
equations (C4) (with the use of complex notation) and (Cll) for P 
using the boundary condition 

P c 

— = Zn at x = 0 (station C) 

Up 1 ° 


gives 


P = 2P r 


Zp cosh(a + jp)x + Z Q sirih(a + j0)x 


u ° Zo + Zc 

where the attenuation constant is the positive root of 




ao 2 p| + R 2 - cop e 


■) 


1/2 


the phase constant is 


P = V- 2 TT + a 2 


K 


(CIO) 


(Cll) 

'distributed 


(C12) 


(C13) 


is, solving 
and U and 


(C14) 


(C15) 


(C16) 


(C17) 
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II 


the specific characteristic impedance of the line is 

Zq = £ (p - ja) (C18) 

CO 

and P q is a constant. (Physically, P q is the disturbance complex pressure 
amplitude at x = 0 for the component of the wave that is traveling in the 
negative x-direction. ) 

Transfer Function 

From equation (C15) the complex amplitude of the pressure perturbation at 
station C (x = 0) is 


P C = ^0 (C19) 

Dividing equation (C15) by equation (C19) gives 

Zs 

~ = cosh(a + jp)x + 7=2. sinh(a + jp)x (C20) 

PC 

from which equations (l) and (2) immediately follow. 


Calculation of Parameters 


In order to calculate the theoretical values of a, p, and Z Q , I x and 
Iy must be found. It Is noted’ that 


N = 



(C21) 


Equation (C2l) and the definition of Bessel-Kelvin functions of the first kind 
given in reference 12 yield 

J 0 (aW) = ber 0 (-a^) + jbei 0 (-a^ j (C22) 

and 

J 1 (aN) = beri(-a ^£p) + jbei-^-aySp j (C23) 

where berQ(-a^po)/p) , beio( -a->/pa>/p) , berq( -a-^/pcD/p) , and beiq( -a-^/pco/jj.) are 
real. Tables (refs. 12 and 13) do not explicitly give values of ber]_ and 
bei^. They can be evaluated from the relations (ref. 13) 

ber-L = _L_ ( ber^ - bei^) (C24) 

VS" 

and 
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i 


beil = ~^z (ber o + bei 0 ) 

Series forms of "berQ, "beiQ, tero, and 13610, given in reference 12, 
tions (C22) to (C25) yield 


(C25) 


and equa- 


Jo(aW) = "bero ^a ^ ^ jtieio ^ ^ 

J x (aN) = -'ber 1 ^a^^j - jbeii 

Combining equations (CIO), (C26), and (C27) and rearranging terms gives 


(C26) 

(C27) 



The values of a, 0, and Z Q can now be obtained for a given frequency from 
tables (refs. 12 and 13) and the appropriate preceding equations. 

Simpler approximate formulas for the attenuation constant a, the phase 
constant 0, and the specific characteristic impedance of the line Zo are 
given in references 7 and 10. 
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APPENDIX D 


PHYSICAL CONSTANTS OF TEST OIL 

The values of the density p of the test oil were obtained from the manu- 
facturer. All other physical constants were obtained from references 14 and 15 
for an oil of the same density and viscosity as those of the test oil. The 
minimum Prandtl number for the oil in the tests was 28; the ratio of specific 
heats in the experiment was 1.159. The density and bulk .mod-ulus of the oil are 
given in the following table : 


Run 

Temperature 

Density 

Adiabatic bulk modulus 8. 


°F 

°K 

slugs /cu ft 

kg/cu m 

lb/sq ft 

N/sq m 

a 

125 

325 

1.65 

850 

3. 04x10 7 

1. 46x10 9 

b 

100 

311 

1.67 

861 

3.22 

1.54 

c 

80 

300 

1.68 

866 

3.37 

1.61 


a Corrected for line compliance as in ref. 1. 
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APPENDIX E 


DETERMINATION OF BEST CURVE THROUGH EXPERIMENTAL 
ATTENUATION CONSTANT DATA POINTS 
A "best-fit curve of the form 

ajp = A -y/f 

is determined by a least-squares method. The constant A is chosen such 


- A y? i J 2 = Minimum 


Equation (E2) requires 


n , 


At/fi jWfi = 0 


Solving equation (E3) for the desired value of A gives 

L a i yU 


A = 


£ *i 


(El) 

that 

(E2) 


(E3) 


(E4) 
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